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ABSTRACT 

We examine local Lagrangian approximations for the gravitational evolution of the den- 
sity distribution function. In these approximations, the final density at a Lagrangian point q 
at a time t is taken to be a function only of t and of the initial density at the same Lagrangian 
point. A general expression is given for the evolved density distribution function for such 
approximations, and we show that the vertex generating function for a local Lagrangian 
mapping applied to an initially Gaussian density field bears a simple relation to the map- 
ping itself. Using this result, we design a local Lagrangian mapping which reproduces nearly 
exactly the hierarchical amplitudes given by perturbation theory for gravitational evolution. 
When extended to smoothed density fields and applied to Gaussian initial conditions, this 
mapping produces a final density distribution function in excellent agreement with full nu- 
merical simulations of gravitational clustering. We also examine the application of these 
local Lagrangian approximations to non-Gaussian initial conditions. 
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1 INTRODUCTION 



A variety of statistics have been developed to describe the evolution of the large-scale 
matter distribution in the universe. Among these is the one-point probability distribution 
function (PDF) of the density field, P(p), which gives the probability that the density at a 
random point in space lies between p and p+dp. In the linear regime, the entire density field is 
simply scaled up the growth factor D(t), so P(p) (with p suitably rescaled) remains constant. 
In the nonlinear regime, however, P[p) evolves in a complex manner. A number of recent 
studies have examined the evolution of the PDF of the density field during gravitational 
clustering in the nonlinear regime. Kofman (1991), Kofman et al. (1994), and Bernardeau 
& Kofman (1995) used the Zel'dovich approximation and the condition of mass conservation 
to derive an approximate expression for P(p). Padmanabhan & Subramanian (1993) derived 
an approximation to the smoothed PDF from the Zel'dovich approximation. Juszkiewicz et 
al. (1994) used the Edgeworth expansion along with the moments of the evolved distribution 
to obtain an approximation for the evolved PDF (see also Bernardeau & Kofman 1995.) [For 
a recent review of approximation methods in general, see Sahni & Coles 1995]. 

All of these approximations assume Gaussian initial conditions; much less is known about 
the evolution of the PDF for non-Gaussian initial conditions. A number of numerical sim- 
ulations have been performed to investigate the evolution of the density field in various 
non-Gaussian toy models (Messina et al. 1990; Moscardini et al. 1991; Matarrese et al. 
1991; Weinberg & Cole 1992; Coles et al. 1993). Fry & Scherrer (1994) examined ana- 
lytically the evolution of skewness in arbitrary non-Gaussian models, and this analysis was 
extended to the kurtosis by Chodorowski & Bouchet (1996), but very little analytic work 
has been done on the evolution of the full PDF in arbitrary non-Gaussian models. 

In this paper, we examine a class of approximations for the evolution of the PDF which 
we call "local Lagrangian approximations" . By "local Lagrangian" , we mean that the density 
at the Lagrangian point q at a time t is approximated as a function only of t and the initial 
value of p(q): 

p(q,*) = /(A>(q),*), (1) 

where po(q) = p(q, to). [We caution the reader that the term "local" has been used in this 
context with a variety of different meanings]. The linear approximation, for example, is a 
local Eulerian mapping. The Zel'dovich approximation is a Lagrangian mapping in which the 
density at q is a function only of t and the second partial derivatives of the initial potential 
00 (q): 

p(q,t) = /(^,t) (2) 

For certain special cases (1-dimensional collapse, spherical collapse), the Zel'dovich approxi- 
mation reduces to a local Lagrangian mapping of the initial density field as given in equation 
(H]). In fact, these cases form the basis of two of our approximations. Note that we do not 
require a prescription for the Lagrangian mapping x(q); all we care about is the function 
given in equation ([!]), which is sufficient to calculate the PDF P(p). In fact, for a given 
mapping / in equation (]!]), there may not even be a Lagrangian mapping scheme which pro- 
duces /; nonetheless our mapping could still provide a good approximation to the evolution 
of the PDF. 

Our motivation for considering such mappings is two-fold: for the case of Gaussian initial 
conditions, it would be extremely interesting if a simple mapping of the form given in equation 
(H|) could give an accurate description of the evolution of the density PDF. We will see that 
this is indeed the case. For non-Gaussian initial conditions, such a mapping would allow for 
the investigation of the evolution of the density PDF for a wide range of initial conditions, 
and it might help to answer some general questions about the evolution of such density 
fields. For example, in the quasi-linear regime, is there any general difference between the 
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rate at which (5 2 ) evolves in non-Gaussian models, versus its evolution in Gaussian models? 
Fry & Scherrer (1994) suggested that in the quasi-linear regime, models with positive initial 
skewness would have an rms density fluctuation slightly larger than that predicted by linear 
theory, while negative skewness models would have an rms fluctuation slightly smaller than 
the linear prediction. The results of Weinberg & Cole (1992) support this conclusion with 
regard to models with negative initial skewness but are inconclusive with regard to models 
with positive initial skewness. A second issue addressed by Fry & Scherrer was the evolution 
of skewness in non-Gaussian models. Their results suggest that the skewness of the evolved 
density field is sensitive to the initial kurtosis (as well as the initial skewness), but the 
expressions they derive contain integrals over the initial three- and four-point functions, so 
it is difficult to draw general conclusions about the evolution of the skewness in such models. 
These are some of the many important general questions about the evolution of non-Gaussian 
models which are at present unanswered. 

In this paper, we consider three different models for the Lagrangian density mapping in 
eq. ([!]). The first of these corresponds to the exact evolution in the 1-dimensional case, 
and it would be exact in the case of 1-D symmetry. Our second model is based on the 
approximation method given by Padmanabhan & Subramanian (1993), and it corresponds 
to the Zel'dovich approximation for spherically-symmetric collapse. We believe that the 
correct PDF will lie somewhere in between these two cases. The third approximation, in 
fact, does lie between these two extremes and is constructed to give approximately the correct 
hierarchical amplitudes in the quasi-linear regime. 

In the next section, we present the motivation for our three approximations and derive 
the form of the density PDF for local Lagrangian mappings. In Section 3, we use our approx- 
imations to derive analytic perturbative results in the quasi-linear regime for the evolution 
of both Gaussian and non-Gaussian initial conditions, and we discuss the effects of smooth- 
ing. In Section 4, we use our approximations in the nonlinear regime to derive the evolved 
PDF for Gaussian initial conditions, which we compare with numerical simulations. We then 
apply our approximations to a variety of non-Gaussian initial conditions. Our results and 
conclusions are summarized in Section 5. In the Appendix, we derive the relation between 
any local Lagrangian mapping of an initially Gaussian density field and the corresponding 
vertex generating function. 



2 LOCAL LAGRANGIAN APPROXIMATIONS 

To obtain the form of the evolved density PDF one has to use some approximation scheme 
describing the particle dynamics, as well as to adopt a set of initial conditions on the density 
field. One of the most efficient approximation schemes is the Zel'dovich approximation (ZA) 
(Zel'dovich 1970) which can be thought of as an operator acting on the initial (Lagrangian) 
comoving position q of a particle and yielding its final (Eulerian) comoving position x. 
Specifically, one has: 

x(q,t)=q + D(t)*(q), (3) 

where D(t) is a universal time dependent function proportional to the expansion factor 
a(t) for a flat, pressureless universe, and SE'(q) is proportional to the gradient of the initial 
gravitational potential. 

The evolution of the PDF of the density field in the Zel'dovich approximation has been 
discussed in detail by Kofman et al. (1994); here we briefly summarize their results relevant to 
our work. Consider a particle of mass dm uniformly spread at time to inside an infinitesimal 
Lagrangian volume d q. At some later time t the respective volume has evolved to the 
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Eulerian infinitesimal one d 3 x q . Mass conservation implies: 



p(q,t)d 3 x q = p(q,t )d 3 q, (4) 

In the quasi-linear regime prior to shell-crossing, multi-streaming may be neglected; here we 
make the assumption that multi-streaming is unimportant. 
From equation ([|) one obtains: 

Pi^t) = p\\^\\-\ (5) 
According to equation (|3|) one has then: 

II Qqll 11-^ + D ^ || . (6) 
Taking r\ = p/p, equation (2.3) yields: 

" = \\1-D%W (7) 

If d^i/dqj has eigenvalues Ai, A2, and A3, then the expression for the evolved density in the 
Zel'dovich approximation can be written as: 

/ +\ ^o(q) /oN 



[l-D(t)X 1 ][l-D(t)X 2 ][l-D(t)X s 



The distribution of A can be calculated exactly in the case of Gaussian initial conditions, and 
the results applied to determine the exact form for P(p) (Kofmanetal. 1994). Unfortunately, 
the exact distribution of the eigenvalues of d^i/dqj is not easy to derive for most non- 
Gaussian models, so we consider several possible ways to simplify eq. (f|). 

We look for approximations to equation (§) which are "local", i.e., the right hand side 
is a function only of D(t) and do- To derive approximations of this sort, we note that the 
relation between <5 and the X's is 

So = Ai + A 2 + A 3 (9) 

Consider first the approximation: 

" (q - i)= [l-D( f )(A 1 + A 2 + A 3 )r (10) 
With equation (§), this reduces to 

[1 - D{t)8o[q)\ 

This approximation is exact in the limit of 1-dimensional collapse (A2 = A3 = 0) and cor- 
responds to exact gravitational evolution in one dimension (Shandarin & Zel'dovich 1989). 
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Hence the results we derive for this approximation will give the exact evolution of the PDF 
in one dimension. Equation (O) has also been investigated as an approximation to the evo- 
lution of the density field in three dimensions (Nusser et al. 1991; Mataresse et al. 1992). 
We will refer to equation ( |Tl~|) as the planar approximation. 

The other extreme case is spherical collapse, for which Ai = A2 = A3. In this case, 
equation (|8|) becomes 

/ j\ 7/0 (q) /, x 

^ = w^wmjw- (12) 

Equation (|I~2"D corresponds to the approximation of Padmanabhan and Subramanian (1993) 
in the limit of zero smoothing; this approximation was also used by Betancort-Rijo (1991) in 
a study of the evolution of the rms density fluctuation. We shall refer to this as the spherical 
approximation. 

We expect that the actual evolution of the density field lies somewhere in between spher- 
ical and planar collapse. Note that both equations fll!]) and (^) are of the form 

l A ^o(q) ( ,o\ 

" (<1 ' t)= (l-D( t )Mq)/«)- <13) 

with a = 1 for the planar approximation and a = 3 for the spherical approximation. For 
our third approximation, we choose a = 3/2 to obtain: 

nfa t) = (14) 

m ' t} [l-2D(t)5 (q)/3]3/2- li4J 

The choice of a = 3/2 has no particular physical significance, but we show in the next 
section [based on earlier results of Bernardeau (1992) and Bernardeau & Kofman (1995)] 
that the hierarichal amplitudes for this model closely mimic the results of exact perturbation 
theory. For that reason, we will refer to equation ([L4]) by the oxymoronic name of the exact 
approximation. 

How do we get from the Lagrangian mappings given by equations ([II]) , fll2|) , and (|14|) to 
the evolved density distribution function? Following Kofman et al. (1994), we define P(p) 
(or, equivalently, P(t])) to be the Eulerian density distribution function and take Q(p) to be 
the Lagrangian distribution function. Basically, P(p) gives the probability that a randomly- 
selected point in space has a density in the interval p to p + dp, while Q(p) is the probability 
that a randomly-selected mass point has a density in that interval. Since the probability 
of "randomly" selecting a given mass point at an Eulerian location x is proportional to the 
density at x, we have (Kofman et al. 1994) Q(p) = (p/p)P(p), or 

P(V) = fav). (15) 

Furthermore, we assume an initial Eulerian distribution Po(^o) an d initial Lagrangian dis- 
tribution Qo(rj ), also related by -Po(?7o) = Qo(f]o)/Vo- Given a Lagrangian mapping of the 
form t) = f(rj (q),t), the evolved Lagrangian PDF is 

Q(v) = Qo(r\v))^-^i (16) 
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where f~ l is the inverse of the Lagrangian mapping given in equation (Jl|). Using the rela- 
tions between the Lagrangian and Eulerian distributions functions, we can express the final 
Eulerian density distribution in terms of the initial Eulerian distribution: 



,-U^r\ri)dr\rj) 



P(v) = Po(f- L (v) r f " ^ • (17) 

All of these equations can be simplified by taking our initial epoch sufficiently early that 
770 ~ 1 and Po(i]o) ~ Qo(Vo)- I n particular, equation (|TTD becomes 

P(v)=Po(r\v)) 1 -^ ]1 . (18) 
r] arj 

Note that in order for P and Q to represent probability distribution functions, they must 
satisfy 

J Q( V )d V = 1, (19) 

and 



J P{ri)dr) = 1. (20) 



Equation ( fl9|) is equivalent to mass conservation and is automatically satisfied by local La- 
grangian mappings of the form given in equation (|l|). Equation (^TJ) gives the conservation 
of Eulerian volume, and it is not automatically satisfied. Assuming that our initial distri- 
bution functions Pq and Qq are correctly normalized, then equation (|20|) is satisfied by the 
Zel'dovich approximation and our planar approximation; it is not satisfied by equation ([13]) 
for a ^ 1, and therefore the spherical and exact approximations fail this test. Hence, we 
cannot use equations ([12]) or (|H| ) as written as valid approximations for the evolution of the 
density. To correct this problem, we modify equation (|]) to read 

r ) (q,t) = N(t)f( m (q),t), (21) 



where N(t) is a time- dependent function given by substituting equation (pT8|) into our nor- 
malization condition [equation (p0|)1: 

m ^nbi PM ^jh) ] (22) 

With the correct normalization, our general mapping given by eq. (|T3|) becomes 

((1 - .D(t)^(q)/ a )°> 
' ( * t)= (l- D ( f )ft(q)/a)«- (23) 

Note that our N(t) resembles the multistreaming factor N s discussed in Kofman et al. (1994), 
although our factor has been introduced as a mathematical construct in order to keep our 
probabilities normalized. 



3 PERTURB ATIVE RESULTS 
3.1 Gaussian Initial Conditions 
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In order to test our mappings for agreement with the true gravitational evolution of 
the PDF, we will first consider the quasilinear case \S(t)\ < 1. In this limit a number of 
perturbative results are known for the evolution of Gaussian initial conditions, and a few for 
non-Gaussian initial conditions, which can be compared with the perturbative predictions 
of our local Lagrangian approximations. 

The density field can be described in terms of the cumulants k p , which are functions of 
the moments of the density field (S p ). The first few cumulants are given by (see Stuart & 
Ord 1987 for a more detailed discussion) 

K2 = (S 2 ) 

% = (s 3 ) 

K4 = (5 4 ) - 3(5 2 ) 2 (24) 

For Gaussian initial conditions, it is possible to show that the cumulants k p of the evolved 

density field satisfy Kp/a 2 ^ 1 ^ — > constant in the limit a — > 0, where a = (S 2 ) 1 ^ 2 , and the 
constants are called the hierarchical amplitudes, denoted S p : 

S p = K> 2(P_1) (25) 

This result was first derived for the skewness (k^) by Peebles (1980), for the kurtosis (K4) 
by Fry (1984), and a method for calculating the full hierarchy of the S p was derived by 
Bernardeau (1992). Our local Lagrangian approximations also produce a hierarchical clus- 
tering pattern (i.e., the cumulants of the evolved distribution satisfy equation ^). Of course, 
there is little point in using approximations to derive S p in the Gaussian case, since these 
values can be calculated exactly. However, it is precisely because the S p s are known exactly 
that the calculation of S p for the Gaussian case can be used to estimate the accuracy and 
general behavior of our approximations. This calculation also serves as a warm-up for the 
case of non-Gaussian initial conditions, for which no general results for S p have been derived 
(although see Fry & Scherrer 1994; Chodorowski & Bouchet 1996). 



First note that all of our mappings (equations |TT|, 12 , [TJ]) can be expressed in the form 



q(q,t) = /ft(q)]> (26) 

where 5i is the linearly-evolved density: Si(q) = D(t)So(q). Consider first the planar mapping 
given by equation flllD . Taking 770 (q) = 1 we obtain: 

tj(q, t) = 1 + Mq) + <Mq) 2 + Mq) 3 + - (27) 

To calculate S3, we need to find k 3 = (S 3 )e, where we now distinguish between Eulerian 
averages (effectively taken over volume and denoted with a subscript E) and Lagrangian 
averages (taken over mass and denoted with a subscript L). Our task is to express the 
Eulerian average of powers of 6 [which are the numbers which enter into equation (|24])1 i n 
terms of the Lagrangian average of powers of 77 [which can be derived from equation (|27|)1. 
For powers of 77, the relation between the Eulerian and Lagrangian averages takes the simple 
form 

(v u )e = (V^L (28) 
which follows directly from equation (0). Then we can express K3 as 



*3 = ((r/-l) 3 ) 



E 



(v 2 )L-m L + 2 (29) 
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Combining equation ( p9] ) with equation (p7|), we obtain, for the planar approximation, 

oo 

«s = £<(;- 2)*?). (30) 

For Gaussian initial conditions, 

$) = U l)!!rr' (jeven), 

= (jodd), (31) 

where a { = (5f) 1/2 = D(t)a is the linearly-evolved rms fluctuation. Then we end up with 

oo 

« 3 = £2(n-l)(2n-l)!!of\ (32) 

n=2 

so that 

oo 

SM) = E 2 ( n ~ !)(2n - l)!!af n - 4 = 6 + 60a? + 0(4), (33) 

n=2 

and S^O) = 6. This result has been derived previously by Bernardeau & Kofman (1995) 
using more complex techniques, but our a 2 term differs from theirs. At this point, we must 
be careful about our definition of the hierarchical amplitudes. In equation ([25]), we have used 
the linearly- evolved rms fluctuation 07, rather than the true rms fluctuation a = (5 2 ) 1/2 . If 
instead we use the true rms fluctuation, we can define a new set of hierarchical amplitudes, 
given by: 

S p (a) = Kp/^y- 1 (34) 

This distinction is usually ignored, because S p (0) = S p (0) for Gaussian initial conditions. 
However, when expanding S p to higher order, or when dealing with non-Gaussian initial 

conditions (Fry & Scherrer 1994), the distinction must be made. Note that S p is closer to 

what observers actually measure. When we expand (S 2 )e — (t))l — 1 5 we obtain 

a 2 = af + 3af + 0(af). (35) 

Substituting this into equation ( |34"|) and reexpressing everything in terms of a rather than 
o~i, we obtain an expression for S p which agrees with Bernardeau & Kofman: 

S 3 (a) = 6 + 24a 2 + 0(a 4 ) (36) 

Note that because o 2 = of + 0(af), it does not matter if we use o~i or a on the right-hand 
side of equation (|36|), although it would make a difference if we expanded out to fourth order 
in a. 

We can similarly expand the expression for K4 to obtain 

«4 = (v 3 )l - 4(v 2 )l - Hv)l + 12(^)l - 6 (37) 

Again, substituting equation ( p7|) for rj in equation (|37j), expanding out term by term, and 
taking the averages appropriate to the Gaussian initial conditions from equation (|3~TD, we 
obtain S^(0) = 72, in agreement with Bernardeau & Kofman (1995). 
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Now consider the spherical approximation (equation O). This approximation is com- 
plicated by the fact that for Gaussian initial conditions! T\l//(?7o, t)) ^ 1, so we have to 
include the normalizing factor given by equation (p2[). For Gaussian initial conditions, the 

ricaTap 



normalizing factor in equation (|22| ) for the sphericaTapproximation is 

N=((l-D(t)5 /3) 3 ) = l + af/3 (38) 
and the spherical mapping becomes: 

"«^(^W (39) 



Substituting this mapping into equations (|29| ) and (137|) and taking the appropriate averages 



for the Gaussian initial distribution from equation (|31|) we obtain 



S 3 (Q) = 4 (40) 
Si(0) = 272/9 (41) 

Oddly, these are identical to the hierarchical amplitudes obtained for the full Zel'dovich 
approximation (see, for example, Munshi et al. 1994; Bernardeau & Kofman 1995). In 
fact, we can show that the spherical approximation and the Zeldovich approximation have 
identical values for all of the hierarchical amplitudes S p (0). To do this, we introduce the 
vertex generating function for the density field (Bernardeau 1992) 

Gs(r) = £ £t» (42) 

where 



n=l 



f(^")(x)^(x 1 )..jW(x n ))/x rf 3 Xl ..A n 

and 5^ is the n^-order expansion of 5. [Our convention for the sign of r is the same as that 
of Munshi et al. (1994) and the opposite of Bernardeau (1992) and Bernardeau & Kofman 
(1995)]. The values of the v n totally determine S p (0), through the relations: 

S 3 {0) = 3z/ 2 , 

S 4 (0) = 4z/ 3 + 12i/ 2 2 , (44) 

and so on (Bernardeau 1992). 

In the Appendix, we demonstrate an extremely useful result for local Lagrangian map- 
pings with Gaussian initial conditions: for a given mapping r/(q, t) = N(t)f[Si(q)], where 
5i is the linearly-evolved density, the vertex generation function is just given by the same 
Lagrangian mapping without the normalizing function: 

G 5 (t) = f(r) - 1 (45) 

[In fact, the same result holds for local Eulerian mappings, but these are not the subject of 
this paper]. For the spherical approximation, equation fl45|) gives 

ft M = JT^JW - 1 (46) 
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which is identical to Gs{r) for the Zeldovich approximation (Munshi, et al.). Thus, S p (0) for 
our spherical approximation and S p (0) for the Zeldovich approximation will be identical for 
all p. This does not mean that the spherical approximation and the Zeldovich approximation 
produce identical evolved PDF's. The reason is that although S p (0) is the same for these 
two approximations, the S p (a)'s differ at higher order in a. The importance of considering 
higher-order a terms has been emphasized by Bernardeau & Kofman (1995), and Scoccimarro 
& Frieman (1996) have recently calculated the higher-order corrections for the Zeldovich 
approximation. Scoccimarro & Frieman obtain, for £3(0"): 

SM) = 4 + (1112/75)*? + 0(of) (47) 
S 3 (a) = 4 + (352/75)a 2 + 0(a 4 ) (48) 

In comparison, our spherical approximation gives: 

S- 3 (<7,) = 4 + (286/27)of + 0(af) (49) 
S 3 (a) = 4+(118/27)o 2 + 0(a 4 ). (50) 

Although the higher-order contributions to S3 are of similar magnitude in the two cases, 
they are not identically equal, so the spherical approximation and Zeldovich approximation 
give different density distributions. This result demonstrates the importance of higher-order 
calculations: two density fields can have identical S p (0) for all p and yet have different PDF's. 

The result given in equation ( (45]) leads naturally to our "exact" approximation (equation 
|i~4"P . An expression for Gg(r) can be calculated in parametric form for the case of the exact 
evolution of the density field (Bernardeau 1992; Munshi et al. 1994; Bernardeau & Kofman 
1995): 

9 (6- sing) 2 
Gs{T) = 2 (1 - 008(9)3 - 1 ' 



for r > 0, and 



Gs(r) 



3 r 3 

5I4 1 



9 (sinh 9 



sm8)} 2/3 , (51) 



2 (cosh g- 1) ; 



T 



- 3 -[ 3 -(smh6-6)f/\ (52) 



for r < 0. Bernardeau (1992) and Bernardeau & Kofman (1995) have noted that an excellent 
approximation for G${t) for the case of exact evolution is given by 

G 5 (r) = (l-2r/3)- 3 / 2 -l. (53) 
Now we can argue backwards from our result in equation (f4~5D: the local Lagrangian mapping 
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will produce a hierarchy of S p (0) very close to the results of exact evolution. Hence, the local 
Lagrangian mapping given in equation ([54]) should provide a good approximation to the true 
evolution in the quasi-linear regime. As expected, the lowest order hierarchical amplitudes 
for this approximation are in excellent agreement with the exact values: our approximation 
gives S 3 (0) = 5, £4(0) = 440/9 « 48.9, compared to S 3 (0) = 4.9 and 5 4 (0) = 45.9 for the 
exact quasi-linear evolution. 

We note in passing that using equations (|5T|)-(|5"2|), along with equation (|45|), it is possible 



to derive a local Lagrangian approximation in parametric form which exactlyreproduces the 
hierarchical amplitudes S p (0) for all p. However, the extra complexity involved in the para- 
metric representation, plus the fact that higher order terms will diverge from their correct 
values anyway, probably makes this approximation less useful than our "exact" approxima- 
tion. 

Now let us evaluate the usefulness of the approximations which we have derived. A 
comparison of Gs(r) for the planar and spherical approximations with G$(t) for the case 
of exact evolution, expanded out in a power series (Munshi et al. 1994) shows that the 
coefficients in the expansion for the planar case are larger than the corresponding coefficients 
for exact evolution, while the opposite is true for the spherical approximation. Thus, S p (0) 
for the planar approximation gives an upper bound on the true S p (0), while 5*^(0) for the 
spherical approximation gives a lower bound. In that sense, these two approximations bound 
the true evolution of the PDF: the planar approximation gives a PDF which deviates more 
strongly from a Gaussian than the true evolution, while the spherical approximation yields 
a PDF which deviates less from a Gaussian. This makes physical sense, since these two 
approximations correspond to planar collapse and spherical collapse, respectively, and the 
true evolution should lie somewhere in between. The "exact" approximation, on the other 
hand, should provide a reasonable approximation to the evolution of the PDF for quasi-linear 
evolution. 

All of these calculations are valid only for Gaussian initial conditions. However, it is 
plausible that our three approximations can be extended to provide some insight into the 
evolution of the PDF for the case of non-Gaussian initial conditions: arguing in analogy 
with the Gaussian case, we expect the planar and spherical approximations to provide up- 
per and lower bounds on the deviation of the PDF from the initial conditions, while the 
"exact" approximation should give a good estimate of the overall evolution. More impor- 
tantly, any features in the evolution shared by all three approximations are likely to be true 
characteristics of the gravitational evolution of the PDF. 

3.2 Effects of Smoothing 

Before venturing into the murky world of non-Gaussian initial conditions, we must con- 
sider the effects of smoothing. Our approximations and results in the previous section apply 
only to the unsmoothed density field, while it is the smoothed density field which is actually 
observed. Since our local Lagrangian approximations give only the PDF and do not provide 
a prescription for actually moving the matter around, there is in principle no way to derive 
smoothed versions of them. However, we can derive plausible "smoothed" approximations 
which give the same hierarchical amplitudes as the smoothed mappings would. 

Our argument is based on the calculations of Bernardeau (1994), who showed that for 
a spherical tophat window function, there is a simple relation between G$(t) for a partic- 
ular density field, and Gf(r) for the corresponding smoothed density field. For simplicity, 
we consider only a density field with a power-law power spectrum P(k) oc k n , for which 
Bernardeau (1994) obtained the implicit equation 

Gf(r) = G 5 (r[l + Gf(r)]-(« +3 )/ 6 ). (55) 
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n 



Consider the vertex function of the density field produced by first applying a given local La 
grangian mapping, and then smoothing. Using equation (p5|), we can generate a "smoothed 
local Lagrangian mapping which produces a final density field with the same vertex generat- 
ing function. To do this, we simply use equation (f45|): for local Lagrangian approximations, 
a transformation of the vertex generating function corresponds to the same transformation 
of the Lagrangian mapping. Thus, for a given mapping r/(q, t) = N(t)f(Si), the "smoothed" 
mapping is given by 

Vs = f(Si) = f(5i[fs(Si)}- {n+3)/6 ) (56) 



where f (Si) must then be multiplied by the normalizing factor specified by equation (pOQ. 
Note that this does not mean that f s (Si) corresponds to the density field derived by applying 
the mapping 77 = f(6i) and then smoothing; rather, fs(Si) will produce a density field with 
the same hierarchical amplitudes S p (0) as the field derived by first applying f(Si) and then 
smoothing. Hence, we expect f (Si) given by equation fl5"6|) to provide a good approximation 
to the PDF of the smoothed density field. 

We will confine our attention to the case n — — 1, because it corresponds roughly to the 
the slope of the CDM power spectrum on galaxy clustering scales, and this choice allows 
us to compare with previous work. For the spherical and "exact" approximations, this case 
produces particularly simple smoothed mappings. For the spherical mapping, equation ( |56"1) 
applied to the mapping given by equation (0) gives (Bernardeau & Kofman 1995) 

rj s = (l + V3) 3 ((l + <V3)- 3 > (57) 
while for the exact approximation, equations ( p^ ) and flT4]) yield (Bernardeau 1994): 



Vs = (6^3 + ^l + 5f/9) 3 ((5i/3 + y/l + 5f/9)- 3 ) (58) 
For the planar approximation, we obtain a cubic equation for rjs, which yields 

Vs = [9 1/3 + Sf/gg 1 / 3 + 5i/3f([g^ + Sf/dg 1 / 3 + Si/3}~ 3 ), (59) 
where g(5i) is defined by 

(/ = l/2 + 5f /27 + (125f + 81) 1/2 /18. (60) 

We can now apply both the smoothed and unsmoothed approximations to the case of non- 
Gaussian initial conditions. 

3.3 Non-Gaussian Initial Conditions 

We now repeat our perturbative calculations from Section 3.1 for the case of non-Gaussian 
initial conditions. The methods are identical to those used in Section 3.1; the only difference 
is that for non-Gaussian initial conditions, all terms of the form (Sf) must be retained, rather 
than being reduced to various powers of a as we did for Gaussian initial conditions (equation 



31]) . This alters not just our final expression derived from our mappings, but also changes 
the normalization factor for each mapping. 

Consider first the rms fluctuation (S 2 ). The most interesting question one can ask about 
(5 2 ) is whether the first correction gives growth which is faster or slower than linear. For 
Gaussian initial conditions, the first correction to linear theory is of order a 4 . This correction 
has recently been examined in detail by Lokas et al. (1996) and Scoccimarro & Frieman 
(1996); they find that for unsmoothed density fields, 

(5 2 ) = of + 1.8af (61) 
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For smoothed initial conditions, the sign of the <r 4 term depends on the initial power index; 
it is negative for n > — 1 but positive for n = —2. 

As noted by Fry & Scherrer (1994), the first correction to (S 2 ) for non-Gaussian initial 
conditions is of order a 3 , rather than a 4 , suggesting that non-Gaussian initial conditions 
should give an earlier divergence from linear behavior. Fry & Scherrer (1994) found that for 
arbitrary initial conditions, 

(5 2 )=v? + f 1 (5?) + pU (62) 

where I [Co] is an integral over the initial three-point function of the non-Gaussian distribu- 
tion. If the second term in equation (|62"D dominates the third term or has the same sign 
as the third term, then the sign of the initial skewness determines whether the rms density 
evolves more or less rapidly than linear, with positive skewness models evolving more rapidly 
and negative skewness models less rapidly. However, all of these results apply only to the 
unsmoothed density field. 

Using our approximations from Section 2 and our perturbative methods from Section 3, 
we obtain, for the unsmoothed density field, the following expressions for (5 2 ): 

(5 2 ) = af + (5f) + (5f)+0(a 5 ), (planar), (63) 
9 53 Fi 

(s 2 ) = ^ + s {6f) + W 6f) + ^ + 0{a5) ' (exact) ' (64) 

(5 2 ) = al + ±(5?) + ^(5t) + ~at + 0(a 5 ), (spherical). (65) 

For the Gaussian case, we see that the first correction term is 3<r 4 for the planar case, 1.61<r 4 
for the exact approximation, and 0.78cx 4 for the spherical approximation. Once more, we see 
that the planar and spherical cases bracket the exact perturbative result, while the exact 
approximation comes close to the exact perturbative value. For the non-Gaussian case, the 
expression for the exact evolution is non-local, as shown in equation (|62|). However, in 
the limit where the long-range correlations in the initial density field are small, equation 
flS2|) reduces to (5 2 ) = (Sf) + (13/21) (Sf). The nonlinear correction term in this case is 
nearly identical to the correction term in the exact approximation, (2/3) (Sf), and is again 
bracketed by the correction terms for the planar and spherical approximations. In all three of 
our approximations, the lowest-order nonlinear correction is a positive multiple of the initial 
skewness, so that positive-skewness models have larger rms densities and negative-skewness 
models have smaller rms densities than predicted by linear theory, in agreement with the 
conjecture of Fry & Scherrer (1994). 

The simulations of Weinberg and Cole (1992) do not unambiguously support this con- 
clusion for the case of smoothed density fields, so we consider what happens when we use 



our smoothed mappings from the previous section. Using the mappings for n = — 1 derived 
in the previous section, we find: 

(5 2 ) = of + lffl + lffl + lot + Oio*), (planar), (66) 
(5 2 ) = af + ^(6t) + y + 0(a 5 ), (exact), (67) 
(S 2 ) = a 2 - l -(5 3 ) + ^(5t) + 2 -at + 0{^), (spherical). (68) 
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We see that the effect of smoothing for n = —1 is to significantly reduce the growth of (5 2 ) as 
compared with the unsmoothed case. In particular, we can draw no unambiguous conclusions 
regarding the dependence of (5 2 ) on the sign of the skewness; our three approximations give 
different answers, and in our exact approximation, there is no skewness dependence at all. 

Now consider the evolution of the skewness itself. As noted in section 3.1, we must distin- 
guish between S 3 (a) = /t 3 /(5f) 4 and S 3 ((t) = /t 3 /(<5 2 ) 4 . For non-Gaussian initial conditions 
these two quantities are not equal even in the limit a — > 0, because of the terms containing 
(Sf) in the expansion for (S 2 ). Here we will consider only S3 (0), because this is the quantity 
measured by observers; it is also the definition of skewness used by Fry & Scherrer (1994). 

Again, using the mappings in Sections 2 and 3.2 along with the expression for S3 in section 
3.1, we can obtain expressions for S^O) for our various approximations. For comparison 
with the Gaussian case (and with previous work) it is convenient to express these results 
in terms of the linearly-evolved kurtosis, defined by K41 =— (Sf) — 3<rf, which vanishes for 
Gaussian initial conditions. For the unsmoothed case, we obtain: 

(69) 
(70) 
(71) 



S 3 (0) = 


(Sf) 


+ 6 + 


V 


JSf) 2 


(planar) 


S 3 (0) = 


(Sf) 

of 


+ 5 + 


3 K 4l 

2af 


4 (5f) 2 
3 af ' 


(exact) 


S 3 (0) = 


(Sf) 

4 


+ 4 + 


K41 
°t 


2 (Sf ? 

3 of ' 


(spherical) 



The exact perturbative result for S^O) for non-Gaussian initial conditions (Fry & Scherrer 
1994) contains non-local terms involving integrals over the the initial three- and four-point 
functions. However, in the limit where these terms are small (e.g., for the case of weak 
correlations in the initial density field), the exact perturbative result is: 



Sm = <f + « + - H?fl, {mct) . ,72) 



This expression is remarkably close to our result for the "exact approximation", and it is 
bracketed by the results for the planar and spherical approximations. 
For the smoothed case, we obtain the results: 

Ss(0) = ^ + 4 + ^-|M!, (planar) (73) 

S 3 (0) = ^ + 3 + ^, (exact) (74) 

S 3 (0) = ^ + 2 + |^, (spherical). (75) 

There are no exact perturbative results with which we can compare our smoothed expressions 
for S3, but they display qualitatively the correct behavior; the effect of smoothing is to reduce 
the evolved skewness. 

The results of this section for non-Gaussian initial conditions confirm the conclusions 
we reached in Section 3.1. Our "exact" approximation for the unsmoothed case agrees well 
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with perturbative results for non-Gaussian initial conditions, while the planar and spherical 
results bracket the known perturbative results. This suggests that our "exact" approximation 
should give reasonable results when we go on to calculate the full evolved PDF for a variety 
of non-Gaussian initial conditions, while the planar and spherical approximations may serve 
as useful bounds on the evolution. The situation is murkier with regard to our smoothed 
approximations, but they show at least qualitatively the correct behavior. 



4 THE EVOLVED PDF 

We can now apply our local Lagrangian mappings to derive the evolved PDF for any 
set of initial conditions using equation fl!8j)- Since we no longer assume Si <C 1, we have to 
modify our generic mapping given in eq. (P3j). Following Kofman et al. (1994) we take 

idfe (76) 

where the normalization factor N(t) is now given by 

tf(t) = (|l-tya|°> (77) 

and the average is taken over the distribution of Si. 

Consider first the case of Gaussian initial conditions. Then for the unsmoothed case, 
equation fll8f) gives 

P(j])dr] = ( - )~ 1 /"" 2 J e -« 2 [l-(V^r 1/Q ] 2 /2<x? + e -a 2 [l+(r 7 /JV)-V«)]2 /2(T 2j ^ 

where N is given by 

N = r — ^ e -W|l - ^\ a d8 h (79) 
J-oo \ / 2irai ot 

Note that this expression for P{rf) bears some resemblence to the distribution function de- 
rived by Bernardeau (1994), [eq. (19)], using more sophisticated techniques, but the two 
distributions are different. 

For our smoothed mappings, the results are even simpler. If we assume a mapping of the 
form in equation (|76f) smoothed with a spherical tophat window function, with a k n power 
spectrum, we obtain from equation ( [Sif) the following relation between the evolved value for 
7] and the linear perturbation Sf. 



a 



'V_\(n+3)/6 _ /V_\-l/a+{n+3)/6 

' N y N J 



10) 



which can be substituted into the Gaussian expression for Si, with the appropriate normal- 
izing factor, to obtain P(rj). For example, for n = — 1 we obtain 



(81) 

For the "exact" approximation (a = 3/2), this reduces to 
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where N in this case is given by 



N 



-Sf/2af 



2na l 



:Si + \l + 



9 



d5 l 



(83) 



We can compare these results with numerical simulations of gravitational clustering. We 
have used the simulations of Weinberg (Juszkiewicz et al. 1995; Lokas et al. 1995) with 
power-law initial conditions, smoothed with a spherical tophat. In Figure 1, we show the 
case n = —1, for o\ = 0.5, 0.75, and 1. The points are the results of averaging eight numerical 
simulations with l-a error bars, and the solid curve is the smoothed "exact" approximation, 
equation ([82]). For o\ = 0.5, the agreement is remarkable, particularly given the simplicity 
of the approximation which led to equation d82|) . There is still reasonable agreement at 
oi = 0.75, but the approximation begins to break down at this point, and agreement is poor 
for u i = 1.0. These results are not surprising, since the "exact" approximation was designed 
to mimic the exact hierarchical amplitudes in the limit where o\ — > 0, and it can be expected 
to break down when the contributions to S p of higher order in o\ become important. In 
fact, Scoccimarro & Frieman (1996) have argued that this should occur at of ~ 1/2, which 
is precisely when our approximation appears to break down in Fig. 1. In Fig. 2, we show 
the planar and spherical approximations along with the numerical simulations for o\ = 0.5 
and o\ = 0.75. For o\ = 0.5, the planar and spherical approximations do indeed bound the 
numerical results, but the spherical approximation breaks down completely at oi = 0.75. In 
fact, it is easy to show from equation ( |31~D that the spherical approximation breaks down (in 
the sense of no longer having central maximum near r\ = 1) at of = 9/20. Again, this is near 
the point at which Scoccimarro and Frieman (1996) predict a breakdown of perturbation 
theory. 

The importance of using our smoothing method derived from Bernardeau (1994) is il- 
lustrated in Fig. 3. There we show, for o\ = 0.5, all three of our approximations for the 
unsmoothed PDF. All three of them differ strongly from the numerical results, showing 
greater deviation from the original Gaussian. This is not surprising, since these approxima- 
tions are only appropriate for the PDF measured for an unsmoothed density field, or for a 
smoothed density field with initial power spectrum P(k) oc k~ 3 . 

Finally, we note that our results in the Appendix indicate that there is also a local 
Eulerian mapping which has a vertex generating function equal to Gs(t) = 1/(1 — 2r/3) 3 / 2 
and which should therefore also give an excellent approximation to the evolved PDF in the 
quasilinear regime. This Eulerian approximation is given by taking r^(x, t) = f[5i(x)}, with 
the mapping given by equation (^) or (|80) for the unsmoothed and smoothed density fields, 
respectively. The final result is a PDF which looks the same as equation (|78|) or equation 

if 

Eu 



but missing the factor of 1/rj which transforms the Lagrangian PDF Q{rf) into the 
erian PDF P(r]). So, for example, the final P(r)) for the Eulerian version of the "exact" 
approximation smoothed with a tophat window function for n = — 1 is simply obtained by 
multiplying equation (p2[) by r\: 



P(r])dr] 



■Kir 213 



+ 



(lr 4/ vp( 



9 



2V2W N N' 
where the normalization factor is now given by: 



[( 



8a? lK N' 



' N J 



-i/3i2\ dr] 
J > N 



N 



-6f/2af 



2-KGi 



;Si + \ 1 + 



dSi 



(85) 
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rather than the expression in equation (|83|). This Eulerian approximation is compared with 
the numerical results for o\ = 0.5 in Figure 4. The agreement between this approximation 
and the numerical simulations is excellent. It may seem implausible that both the Eulerian 
and Lagrangian local mappings could produce nearly the same final P{i]), since the distri- 
butions given by equations ( |8"2] ) and (E3I) differ by a factor of rj. However, this difference is 
compensated by the different values for N used in the two equations. In fact, our argument 
in the appendix indicates that these two approximations should give equally good agreement 
with the evolved PDF in the limit where 07 <C 1. 

Now consider the evolution of non-Gaussian initial conditions. Our approximation can 
be applied in an elementary way to any initial density distribution, but there are an infinite 
set of distributions to choose from. To explore the differences which the initial skewness and 
kurtosis make in the evolution, we have chosen four representative distributions; one each 
with positive and negative skewness, and two symmetric (zero skewness) distributions with 
positive and negative kurtosis. The distributions we examine below represent some extreme 
cases and are not physically motivated. However, they give a general idea of the effects of 
positive and negative skewness or kurtosis on the evolution of the PDF. 

We will consider only the smoothed exact approximation with 07 < 0.75, which we 
know produces results in good agreement with the true PDF for the Gaussian case, and for 
definiteness we will take n = — 1. For this case, the mapping in equation (BO) becomes 



5i = - 



'N' 



'N' 



(86) 



with N given by: 



N 



-Si + \ 1 + — 
3 V 9 



dS, 



(87) 



where P(Si) is the initial (non- Gaussian) distribution which we are evolving. 

For the case of positive skewness, a simple choice is the gamma distribution with zero 
1 '■ 



mean and variance a 2 - 



P(Si)dSi 



,5 l 



where each value of v defines a different gamma distribution. A mirror image negative 
skewness distribution with zero mean can be obtained by simply changing S\ to — Si in 
equation fl88|) . For definiteness, we take v = 3, and our results for positive and negative initial 
skewness are shown in figures 5 and 6, at 07 = 0.2, 0.5, and 0.75. Both density distributions 
show qualitatively the expected evolution, i.e., the development of increasing skewness and 
a large positive tail. For the positive skewness case (figure 5) this does not represent a major 
change in the shape of the distribution function. The effect is more dramatic in figure 6, 
where the distribution function with negative initial skewness flips into a PDF with positive 
skewness, as expected. Oddly, the function with negative initial skewness develops a larger 
tail at large 77 than does the function with positive initial skewness. 

For symmetric distribution functions we have chosen two "extreme" representative mod- 
els, the bilateral exponential: 



P{Si)dS t 



V2a t 



-V2\ Sl \/a ld6 



19) 
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which has large positive kurtosis, and the uniform distribution 



1 



d8i, \Si\ < V3a h 
(otherwise), 



(90) 
(91) 



2y/3ai 



= 



with large negative kurtosis. Neither of these can be considered a realistic initial distribution, 
but they illustrate the effect of large kurtosis for symmetric initial conditions. The evolved 
PDF's for these two models are given in figure 7 (bilateral distribution) and figure 8 (uniform 
distribution). In both cases, the singularities in the initial density distribution remain in 
the evolved PDF. However, despite the extreme nature of the initial distribution functions, 
in both cases the evolved PDF shows the expected quantitative behavior, with increasing 
skewness and the development of a tail at large 77. 



5 CONCLUSIONS 

Despite the simple-mindedness of local Lagrangian approximations, our "exact" approx- 
imation provides remarkable agreement with numerical simulations of the evolution of the 
density distribution function with Gaussian initial conditions. This agreement can be under- 
stood in terms of the fact that this approximation reproduces nearly exactly the hierarchical 
amplitudes at tree level. The planar and spherical approximations appear to bound the 
evolution of the PDF but are much less useful (unless, of course, one is interested in the 
evolution of one-dimensional density fields, in which case the planar approximation is exact 
for any initial conditions). 

For the case of non-Gaussian initial conditions, we cannot be as confident. Unlike the 
case of Gaussian initial conditions, the evolution of the hierarchical amplitudes is non-local, 
as has been shown by Fry & Scherrer (1994) and Chodorowski & Bouchet (1996). Thus, 
no local approximation can exactly reproduce the hierarchial amplitudes for non-Gaussian 
initial conditions. However, the exact approximation does a reasonable job of reproducing the 
hierarchical amplitudes for limiting cases where a local approximation is valid. Furthermore, 
the application of this approximation to various non-Gaussian initial conditions does show 
reasonable agreement with the expected qualitative behavior. 

The mapping which gives us the exact approximation can also be applied backwards, 
to map the evolved distribution function back onto the initial distribution function. This 
procedure is guaranteed to produce the correct initial distribution function only for the case 
of Gaussian initial conditions, and only for reasonably small 07, but it is obvious from Fig. la 
that it would be highly accurate in this case. This method should, in principle, be capable of 
distinguishing Gaussian from non-Gaussian initial conditions, even if it could not accurately 
give the exact form of the latter. 



We thank Josh Frieman for helpful discussions. We are grateful to David Weinberg for 
providing us with results of his gravitational clustering simulations. R.J.S. was supported 
in part by the Department of Energy (DE-AC02-76ER01545). Z.A.M.P. and R.J.S were 
supported in part by NASA (NAG 5-2864). 



APPENDIX: The Vertex Generating Function for Local Approximations 
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Here we demonstrate the simple form of the vertex generating function for local La- 
grangian approximations applied to Gaussian initial conditions; namely, for a local La- 
grangian mapping 7] = N(t)f(Si), where N(t) is a function only of 07, the vertex generating 
function is given by 

G,(r) = /(<*,)- 1 (92) 

Our argument proceeds in two stages. First we show that this relation holds for the simpler 
case of local Eulerian mappings. Then we show that the vertex generating function is the 
same for local Eulerian and local Lagrangian mappings. 
Consider first the local Eulerian approximation 

tj(x) = N E (t)f(S^)). (93) 

The mapping / can be expanded in a power series 

00 

/W=E¥? (94) 
3=0 

where 60 — 1 5 and N(t), for Gaussian initial conditions, can be a function only of af. 

00 

N E (t) = J2c k af. (95) 

k=0 

where cq = 1. These three equations give us the expression for the n th order expansion of 5: 

S (n) = E b jCk 6fa k . (96) 

j-\-k=n 

We can substitute this expression into equation ( fi3|) to obtain u n . When we do this and take 
the connected average, all of the terms vanish except for the term j = n, k = 0; the other 
terms are not connected. The final result is 

/(M( 1 )"(x)^ 1 )(x 1 )...^ 1 )(x n )) c d 3 x rf 3 Xl ..A n 
Vn (J(5( 1 )(x)5( 1 )(x'))d 3 xd 3 x') n 

= b n n\ (97) 

since there are n\ different connected graphs. Then from equation (H2J) we obtain 



g 5 (t) = J2b n 



T 



n=l 



= /(r)-l (98) 

Thus, equation fl^o] ) holds for local Eulerian mappings; we now show that G${r) is the 
same if we take f(Si) to be a local Lagrangian mapping instead of a local Eulerian mapping. 
To do this, we first note that G$ is completely determined by the values of v n (equation , 
which are, in turn, completely determined by the S p (OYs (equation il]) . Thus, it suffices 
to show that the local Eulerian mapping given by equation (|93| ) and the local Lagrangian 
mapping 

rj(q) = iV L (t)/(^(q)) (99) 
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have the same values for S p (0), which is equivalent to the statement that they have identical 
cumulants k p in the limit a — > 0. 

To demonstrate this, we introduce the characteristic function 0s (i) for the Eulerian 
mapping, which is the Fourier transform of the PDF: 

Mt) = J PE(v)e ivt dr) (100) 

where we assume that Pe{v) is the PDF for the Eulerian mapping given by equation fl93|). 
If instead we use the same function to produce the local Lagrangian mapping givenby 
equation (|99|), we obtain the PDF Pl(^), with corresponding characteristic function 4>L,(t)- 
Using equation ([18]) for the Lagrangian mapping, and its equivalent (without the 1/rj factor) 
for the Eulerian mapping, we find that (pi and 4>e are related by: 

M^t) = U' L {t\ (ioi) 

N E i 

where Nl and Ne are the normalizing factors for the Lagrangian and Eulerian mappings, 
given respectively by 

N L = < JL> (102) 

and 

N. = ^ (103) 

The cumulants we wish to calculate are related to the characteristic function by 

00 (it) p 

<P(t) = exp(£ ( —tK p ) (104) 
P =i P ] 

If we let K p E > represent the cumulants for the Eulerian mapping, and K p L ^ represent the 
c umu lants for the Lagrangian mapping, then we can substitute equation (|104|) into equation 
( |i0T|) to obtain: 

ti4 E) ((fWf)Y - 4 L) ] ij T = ln t x + £ ( 105 ) 

P =i P- p=i P- 

where we have used the fact that k[ l ^ = 1. Now we note that for both the Eulerian and 
Lagrangian mappings considered here, the "evolved" distribution is hierarchical, so that 
4 L) ~ 0(a 2( - p -V), and ^ ~ 0{a 2 ^). Furthermore, {(f)(l/f)) p = l + 0(a 2 ). Expanding 
out the right-hand side of equation (|105|) and equating terms with equal powers of t p , we 
find that 

4 E) - 4 L) ~ 0(a 2p ). (106) 

Furthermore, k p for P{rj) and k v for P(S) are identical for p > 1, since rj and 5 differ by a 
constant. Thus, the difference between S p (a) for the Lagrangian local mapping and S p (a) 
for the local Eulerian mapping vanishes in the limit where a — > 0, so that the two mappings 
have the same Gs{t). 

These results allow us to generate, in a simple way, two different mappings which, when 
applied to Gaussian initial conditions, yield a density field with any desired vertex generating 
function G$(t) or hierarchical amplitudes S p (0). 
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FIGURE CAPTIONS 



Figure 1: Comparison of our exact approximation for the evolution of the smoothed density 
distribution function (solid curve) with the results of a numerical simulation of gravitational 
clustering (points with error bars), for Gaussian initial conditions with power spectrum 
P(k) oc k~ l and spherical top-hat smoothing. The distribution functions are calculated for 
a linearly evolved rms fluctuation of (a) 07 = 0.5, (b) 07 = 0.75, (c) 07 = 1.0. 

Figure 2: As Figure 1, but here the results of the numerical simulation of gravitational 
clustering (points with error bars) are compared with the smoothed density distributions 
given by the planar approximation (dotted curve) and the spherical approximation (dashed 
curve) for (a) 07 = 0.5, (b) 07 = 0.75. 

Figure 3: As Figure 1, but here the results of the numerical simulation of gravitational 
clustering (points with error bars) are compared with the unsmoothed density distributions 
given by the exact approximation (solid curve), the planar approximation (dotted curve), 
and the spherical approximation (dashed curve), for 07 = 0.5. 

Figure 4: As Figure la, using an Eulerian version of the exact approximation (solid curve). 

Figure 5: Evolved smoothed density distribution functions given by the exact approximation 
with P(k) oc fc -1 , for an initial gamma function density distribution with positive skewness. 
The distribution functions are calculated for a linearly evolved rms fluctuation of 07 = 0.2 
(solid curve), 07 = 0.5 (dashed curve), and 07 = 0.75 (dotted curve). 

Figure 6: As Figure 5, for an initial gamma function with negative skewness. 

Figure 7: As Figure 5, for an initial bilateral exponential function. 

Figure 8: As Figure 5, for an initial uniform distribution. 
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Fig. la 




Fig. lb 
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Fig. lc 
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